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Raman study of lattice dynamics in Weyl semimetal TaAs 

H. W. Liu/’QP- Richard, 1,2 ’^] Z. D. Song, 1 L. X. Zhao, 1 Z. Fang, 1,2 G.-F. Chen, 1,2 and H. Ding 1,2 

1 Beijing National Laboratory for Condensed Matter Physics, 
and Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 
2 Collaborative Innovation Center of Quantum Matter, Beijing, China 
(Dated: June 16, 2015) 

We report a polarized Raman study of Weyl semimetal TaAs. We observe all the optical phonons, 
with energies and symmetries consistent with our first-principles calculations. We detect additional 
excitations assigned to multiple-phonon excitations. These excitations are accompanied by broad 
peaks separated by 140 cm” 1 that are also most likely associated with multiple-phonon excitations. 

We also noticed a sizable Bi component for the spectral background, for which the origin remains 
unclear. 

PACS numbers: 78.30.-j, 63.20.-e 


I. INTRODUCTION 

Although predicted long ago as a solution of the Weyl 
equation describing relativistic spin 1/2 fermions [I], 
Weyl fermions have never been observed experimentally 
in high-energy physics. A few years ago the existence 
of Weyl nodes behaving like Weyl fermions has been 
predicted in condensed matter physics 0 0 , in mate¬ 
rials called Weyl semimetals. An isolated Weyl node 
can be assimilated in momentum space to a monopole 
of topological charge defined by its chirality [3], How¬ 
ever, only very recently did the existence of Weyl nodes 
proven experimentally. In particular, the existence of 
Weyl nodes predicted in TaAs [D has been inferred ex¬ 
perimentally from its topological surface state measured 
by ARPES [3] and from its negative magneto-resistance 
induced by the chiral anomaly [6j, and more directly from 
its bulk electronic band structure determined from bulk- 
sensitive ARPES measurements [7]. Besides the interest¬ 
ing physics derived from its topological surface and the 
topological nature of Weyl nodes, TaAs has a great po¬ 
tential for applications. Indeed, not only it is a binary 
compound with a relatively simple crystal structure, the 
Weyl nodes are associated with bulk properties, in con¬ 
trast to topological insulators for which only the surface 
shows a particular interest. Despite a complete descrip¬ 
tion of the electronic structure of this exotic material, 
there is only little known, both theoretically and experi¬ 
mentally, about the lattice dynamics of TaAs. 

In this paper, we report a Raman study of TaAs. We 
observe all the optical modes of TaAs at the Brillouin 
zone (BZ) center (r), which consist in one Ai mode, two 
Bi modes and three E modes. The experimental ener¬ 
gies of these modes, along with their symmetries, are 
consistent with our first-principles calculations. We de¬ 
tect peaks beyond the range of the single-phonon excita¬ 
tions assigned to multiple-phonon excitations, as well as 
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broader features separated by 140 cm -1 that also likely 
correspond to multiple-phonon excitations. Finally, we 
report a sizable Bi component for the spectral back¬ 
ground, which has an undetermined origin. 


II. EXPERIMENT 

The single crystals of TaAs used in this study were 
grown by chemical vapor transport and characterized by 
X-ray diffraction to determine the crystal orientation [6] . 
Freshly prepared platelike samples with typical size of 
0.4 x 0.4 x 0.08 mm 3 were prepared for Raman scattering 
measurements. The measurements were performed with 
the 514.5 nrn and 488.0 nm excitations of an Ar-Kr laser 
focussed on flat sample surface regions with a 100 x ob¬ 
jective mounted in a back-scattering micro-Raman con¬ 
figuration. The power at the sample was smaller than 
0.4 mW. The signal was analyzed by a Horiba Jobin 
Yvon T64000 spectrometer equipped with a nitrogen- 
cooled CCD camera. 


III. CALCULATIONS 

The non-centrosymmetric structure of TaAs is char¬ 
acterized by the space group I41md (Cj 1 , group no. 
109), with 4 atoms in one unit cell. An analysis in 
terms of the irreducible representations of this group 
shows that the vibration modes of this system decompose 

[8] into [Ai+E]+[Ai+2Bi-|-3E], where the first and sec¬ 
ond terms correspond to the acoustic and optic phonon 
modes, respectively. While all the optic modes are Ra¬ 
man (R) active, only the Ai and E modes are infrared 
(IR) active. To determine the phonon vibrational con¬ 
figurations and estimate the phonon frequencies, we per¬ 
formed calculations using the first-principles pseudopo¬ 
tential plane wave method package Quantum Espresso 

[9] . We set al2xl2xl2 monkhorst-pack momentum 
(k) point mesh and a 45 Ry cutoff for the electronic 
wavefunctions. The exchange and correlation functional 
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was treated within the generalized gradient approxima¬ 
tion (GGA) of Perdew-Burke-Ernzerhof [TOj- Coordi¬ 
nates and cell shape of experimental data m have been 
fully relaxed until the forces acting on atoms are all 
smaller than ICC 4 Ry/aB and the pressure is smaller than 
0.2 kbar. Using the information on the ground state, it is 
easy to use the Phonon package, which implements the 
density functional perturbation theory (DFPT) 
to get the phonon frequencies and vibration modes at 
the r point. All the optic vibration modes of TaAs are 
illustrated in Fig. jTj For each mode we indicate the cor¬ 
responding irreducible representation, the experimental 
and calculated phonon frequencies, as well as the optical 
activity. While the E inodes are related to vibrations in 
the ab plane, the Ai and Bi modes are associated with 
vibrations along the c axis. In addition, we also calcu¬ 
lated the phonon dispersion along high-symmetry lines 
in the momentum space. We first applied DFPT calcu¬ 
lations using the Phonon package on a coarse 4x4x4 
monkharst-pack phonon momentum (g) point mesh to 
get their frequencies and modes. The dynamic matrix in 
real space is then obtained by Fourier-transformation of 
these modes. Finally, we transformed this real space dy¬ 
namic matrix back to general q points and diagonalized 
it to get the dispersion. 


IV. RESULTS AND DISCUSSION 


Experimentally, the symmetry of the vibration modes 
can be described by the relevant Raman tensors of the 
C 4 ,, group (note that the A 2 channel is not Raman ac¬ 
tive): 
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FIG. 1. (Color online). Main atomic displacements for the op¬ 
tical modes of TaAs. The displacements of Ta and As atoms 
are indicated by black and red arrows, respectively. We cau¬ 
tion that the displacements of Ta atoms are not represented 
when they are found too small compared to the displacements 
of As atoms. The first line below each configuration of vibra¬ 
tion indicates the mode symmetry in the C 4 „ group notation. 
Numbers in parenthesis are used to specify modes when there 
are more than one mode with the same symmetry The sec¬ 
ond line below each configuration of vibration indicates the 
experimental(calculated) mode energy. The last line below 
each configuration of vibration gives the optical activity, with 
I = infrared active, R = Raman active and I+R = infrared 
and Raman active. 
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In Fig.[2]we show the Raman spectra obtained in differ¬ 
ent ab-plane configurations of incident (e*) and scattered 
(e s ) polarizations. Pure Ai and Bi symmetries are ob¬ 
tained in the x'x' and x'y' configurations of polarizations, 
respectively. As expected from our group theory analy¬ 
sis, only one mode is observed in the x'x' spectrum. Its 
energy, 251.9 cm -1 , is slightly higher than the calculated 
value (243.2 cm -1 ). Similarly, we detected two peaks 
in the x'y' configuration, in agreement with our analy¬ 
sis. The Bi(l) peak is observed at 171.7 cm -1 while the 
Bi(2) peak is detected at 259.2 cm -1 . In agreement with 
their symmetries, both the A 4 and the two Bx modes are 
observed in the xx configuration while none of them ap¬ 
pears in the xy spectrum. 

In order to confirm further the symmetry of the Ra¬ 
man peaks observed, we performed an angular depen¬ 
dence study of the ab spectra. In Figs. [ 2 ]+) and^e) 


we display the intensity plots of the Raman spectra as 
a function of the in-plane angle 6 defined with respect 
to the x axis for e*||e s and e* _L e s , respectively. The 
intensity of the Ai and Bi modes for the corresponding 
plots are given in Figs. [2](d) and[2](f), respectively. As 
expected, the Ai peak is observed only clearly for e*||e s , 
with an intensity insensitive to the angle, as shown in 
Fig. if). In contrast, the Bi modes exhibit a strong an¬ 
gular dependence for both the e*||e s and e* _L e s configu¬ 
rations. For e*||e s , we expect from the Bi Raman tensor 
that the intensity should vary with angle like |ccos( 2 d)| 2 , 
with a four-fold symmetry characterized by a maximum 
at 0 degrees and a node at 45 degrees, which is consis¬ 
tent with our observation, as shown in Fig. id). For 
e* _L e s , this pattern is shifted by 45 degrees as the Bi 
Raman tensor predicts an angular variation of the inten¬ 
sity proportional to |csin( 20 )| 2 , in agreement with our 
observation (see Fig. [2jf)). Although this cannot be seen 
clearly from the individual spectra, we also point out a 
very small leak of the Ai mode into the image plot dis- 
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FIG. 2. (Color online) (a) Room temperature Raman spectra of TaAs for various ab-plane configurations of the incident and 
scattered light polarizations. The inset illustrates the experimental configuration, (b) Same as (a) but for the ac plane. The 
background of the ac spectra has been removed and the spectrum in the zx configuration has been multiplied by 10 for a better 
visualization, (c) Intensity plot of the Raman shift as a function of the in-plane angle defined with respect to the x axis, for 
e l ||e s . (d) Angular dependence of the intensity of the Raman modes observed in (c) with e‘||e s . (e) Same as (c) but for e'lei 
(f) Same as (d) but for the Raman modes observed in (e) with e l ie 1 . (g) Long-range spectra of TaAs recorded with 488.0 nm 
laser excitation in several ab-plane configurations of polarization. The dominant symmetries of some sharp excitations at high 
energy are indicated. The inset shows the ab-plane angular dependence of the Raman intensity recorded with 514.5 nm laser 
excitation under the e l ||e s and e 1 _L e s configurations of polarization. The red horizontal lines indicate the broad humps at 
240 cm -1 , 380 cm -1 and 520 cm -1 . The contrast has been adapted to illustrate a sizable Bi symmetry component for the 
spectral background. 


played in Fig. [2je), obtained with e* J_ e s . We attribute 
this leak to a small misalignment of the relative angle 
between e* and e s away from perfect right angle. 

We now describe Raman measurements of the ac plane. 
For perfect alignment, no Bi phonon should be observed 
while only the Aj phonon should be detected when e*| |e s . 
This is consistent with the zz spectrum displayed in 
Fig- lb), which exhibits the Ai mode at the same energy 
as in the ab plane spectra. Despite a much smaller inten¬ 
sity, we can detect some peaks in the zx configuration. In 
principle, only E modes should be observed in this config¬ 
uration. However, four peaks (confirmed in several mea¬ 
surements) are detected instead of the three predicted. 
By matching the experimental data with the calculated 
values, we can unambiguously identify the E(l) and E(2) 
modes at 125.6 cm -1 and 232.6 cm -1 , respectively. This 
leaves two candidate peaks at 251.7 cm -1 and 260.9 cm -1 
to be assigned to the E(3) phonon, calculated to be lo¬ 
cated at 244.3 cm -1 , as shown in Fig. [lj The most likely 
explanation for the origin of the 251.7 cm -1 is its assign¬ 
ment to the Ai mode, which has the same energy within 
error bars. Although this would be strictly forbidden for 
cross polarizations, we explained above a possible mis¬ 


alignment of the e* and e s polarizations away from per¬ 
fect right angle. Since the Ai is much stronger in the zz 
spectrum as compared to the Raman intensity recorded 
in the zx spectrum, this hypothesis is quite plausible. 
Following this conclusion, we tentatively assign the peak 
at 260.9 cm -1 to the E(3) phonon. Although this en¬ 
ergy is close to that of the Bi(2) mode (259.2 cm -1 ), the 
energy difference is measurable. 

So far we discussed single-phonon excitations, which 
have decent agreement with our calculations. The up¬ 
per energy limit is 260.9 cm -1 . However, many Raman 
excitations are observed beyond that limit, as shown 
in Fig. (2jg). Considering that the phonon modes dis¬ 
perse in the momentum space, and allowing for small 
discrepancies between experiments and calculations, the 
energy range where these extra excitations are detected is 
consistent with double-phonon excitations, as confirmed 
by our calculation of the phonon dispersion along high- 
symmetry lines displayed in Fig. [3j It is worth noticing 
at this point that these additional excitations are ob¬ 
served for all the samples that we measured, and also 
with the two laser wavelengths excitations used in our 
study. For example, Fig. |2jg) clearly shows their exis- 
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FIG. 3. (Color online). Calculation of the phonon band 
dispersions in TaAs along high-symmetry lines. The inset 
gives the definition of high-symmetry points in the momen¬ 
tum space. The blue arrows show the momentum path corre¬ 
sponding to the dispersion displayed. 


fence under 488.0 nm excitation while its inset shows 
clearly some intensity in the same energy range, as mea¬ 
sured under 514.5 nm excitation. By comparing the spec¬ 
tra in Figs. Ha) and gg), we report a variation of the 
relative intensity of the Ai and Bi phonons as the laser 
excitation varies. In contrast to the measurements per¬ 
formed using a 514.5 nm laser wavelength, for which the 
spectra show a stronger intensity of the Ai peaks as com¬ 
pared to the Bi peaks, the intensity of the Bi peaks be¬ 
comes larger than that of the Aj mode under 488.0 nm 
laser excitation. 

Double-phonon excitations can occur not only at the 
r point, but also anywhere on the BZ boundary, with 
the only restriction that the total momentum is 0. A 
complete assignment of each feature is rather compli¬ 
cated and beyond the scope of our study, not only be¬ 
cause modes disperse but also because their symmetries 
vary with momentum. Nevertheless, the symmetries of 
many excitations can be described in terms of the C± v 
point group and at least two of them can be tentatively 
identified as double excitations involving phonons at F. 
Indeed, we observe two sharp peaks at 358.8 cm -1 and 

386.7 cm -1 that appear clearly in the x'x' and xy spectra. 
Such behavior is expected for a B 2 excitation. We also 
notice some reminiscence of that peak in the x'y' spec¬ 
trum, corresponding to a Ai symmetry. As the ground 
state coincides with the Ai representation, the symme¬ 
try of the excitation is directly determined by the tensor 
product of the irreproducible representations character¬ 


izing the two phonon modes involved. The only three 
possibilities to get a double-phonon excitation at T with 
a B 2 symmetry are: 1) Ai x B 2 = B 2 ; 2) A 2 x B x = B 2 ; 
3) E x E = Ai + A 2 + Bj + B 2 . Since there is no A 2 
and B 2 mode at F, the only possibility is to combine E 
inodes. Interestingly, the excitations at 358.8 cm -1 and 

386.7 cm -1 coincide almost perfectly with [E(l) + E(2)] 
(125.6 cm -1 4 - 232.6 cm -1 = 358.2 cm -1 ) and 
[E(l) + E(3)] (125.6 cm ' 1 + 260.9 cm ' 1 = 386.5 cm” 1 ), 
respectively. 

We also observe additional features most likely related 
to double-phonon excitations, although a precise assign¬ 
ment of the modes involved and the corresponding mo¬ 
mentum locations remains uncertain. Among peaks with 
a dominant but not necessarily exclusive Ai symmetry 
(as described within the Ca v point group), we detect 
features at 468.5 cm -1 , 501.7 cm -1 and 523.1 cm -1 . 
Similarly, the peaks at 490.2 cm -1 , 509.2 cm -1 and 

567.7 cm -1 have a dominant Ai + Bi symmetry. Fi¬ 
nally, we notice some broader features with a full width 
of about 100 cm -1 . The centers of these broad features 
are located at 240 cm -1 , 380 cm -1 and 520 cm -1 , as 
shown in Fig. [2jg) and its inset. Interestingly, these 
broad humps are separated by 140 cm -1 . These fea¬ 
tures are detected in all the ab-plane configurations of 
polarization. 

Interestingly, the inset of Fig. |2jg) suggests that the 
spectral background itself is carries a sizable Bi symme¬ 
try component. Indeed, we observe an intensity modu¬ 
lation with angle that is similar to that of a Bx phonon, 
with lower intensity at 45 degrees for e 4 ||e s , and at 0 
and 90 degrees with e* le s . Whether this phenomenon 
origins from lattice vibration or is related the Weyl 
semimetal nature of the electronic structure of TaAs, is 
unclear. 


V. SUMMARY 

In summary, we performed a polarized Raman study of 
Weyl semimetal TaAs. We identified all the optic phonon 
inodes of this material. Their energies and symmetries 
are consistent with our first-principles calculations. We 
observed additional features assigned to double-phonon 
excitations, as well as much larger humps separated by 
140 cm -1 that are also likely to have a multiple-phonon 
origin. Finally, we identified a sizable Bi component for 
the spectral background, which origin is unknown. 
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